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ABSTRACT 

Using the spectroscopic sample of the SDSS DRl we measure how gas was transformed into 
stars as a function of time and stellar mass: the baryonic conversion tree (BCT). There is a 
clear correlation between early star formation activity and present-day stellar mass: the more 
massive galaxies have formed about 80% of their stars at z > 1, while for the less massive 
ones the value is only about 20%. By comparing the BCT to the dark matter merger tree, 
we find indications that star formation efficiency at z > 1 had to be about a factor of two 
higher than today {^^ 10%) in galaxies with present-day stellar mass larger than 2 x lO^^A/©, 
if this early star formation occurred in the main progenitor Therefore, the LCDM paradigm 
can accommodate a large number of red objects. On the other hand, in galaxies with present- 
day stellar mass less than 10^^ M©, efficient star formation seems to have been triggered 
at z ^ 0.2. We show that there is a characteristic mass (M* ^ 10^" M0) for feedback 
efficiency (or lack of star formation). For galaxies with masses lower than this, feedback 
(or star formation suppression) is very efficient while for higher masses it is not. The BCT, 
determined here for the first time, should be an important observable with which to confront 
theoretical models of galaxy formation. 

Key words: galaxies: fundamental parameters, galaxies: statistics, galaxies: stellar content 



1 INTRODUCTION 

In the current paradigm for galaxy formation, galaxies form in cold 
dark matter halos, which evolve from small, primordial, Gaussian 
fluctuations, by gravitational instability. This mechanism fits well 
in the successful LCDM picture which correctly describes the Uni- 
verse from the Cosmic Microwave Background slI z = 1088 t o 
local galaxy clustering jSpereel et al]l2003t iPercival et al1l200lh . 
One of the strong predictions of this galaxy formation paradigm 
is the typical redshift of dark matter halo formation (i.e. virializa- 
tion) as a function of ha lo mass and cosmological parameters (e.g., 
ISheth&Toniiet]<2004 ). 

Given that cosmological param eters have been tightly con- 
strained (e.g., ISpergel et alj yOO^)), we can reconstruct the 
avera ge dark halo formation history as a function of ma ss 
(e.g.. IPress & Schechteil Il9l4) : ISheth & Tormerl Jl999l 12004 ). 

Naively, the dark matter halo collapse, i.e. virialization, should trig- 
ger baryonic gas transformation into stars; in addition, subsequent 
dark matter mergers should produce star formation episodes. We 
will show that this simple model is not in agreement with observa- 
tions. 

All we can observe is the integrated light of galaxies' stel- 
lar population; thus, to compare the theory prediction for the dark 
matter halos with observations, the process of how baryons are 
transformed into stars needs to be simulated, either through semi- 



analytical recipes or by means of hydro-dynamic N-body simula- 
tions. Since no theory of star formation has been yet established, 
we do not have a fundamental theory that allows us to compute 
from basic principles the star formation efficiency. Further, com- 
plications arise from other phenomena, such as feedback by stars 
and AGN that prevent the formation of giant molecular clouds and 
therefore reduce star formation efficiency. Given the complexity of 
the task of learning about galaxy formation from numerical simu- 
lations, here we take the complementary approach of placing new 
observational constraints on the stellar assembly history as a func- 
tion of galaxy mass. 

In this paper we use about 10'' galaxies from the Sloan Digital 
Sky Survey Data Release 1 (SDSS DRl) to determine, for the first 
time, the amount of baryons that have been transformed into stars 
as a function of total stellar mass and time. This allows us to build 
the baryonic conversion tree (BCT), which can then be compared 
with the dark matter merger tree. We present such comparison and 
show how it can be used to compute the star formation efficiency 
and the relative importance of feedback (or lack of star formation). 

Our main findings are: 

(i) There is a clear correlation between total stellar mass of the 
galaxy and the fraction of gas transformed into stars at z 1. The 
larger the mass the larger the fraction of gas transformed into stars 
at high redshift. 
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(ii) If large galaxies need to be formed by z ■-^ 1 in a "mono - 
lithic" fashion, as obser vations sug g est (e.^ .,^owe^^lJ j^92|): 
[Peacock et al." ('1998); Lillv et al. (1998); Brinchmann 8l EllidI 
(2000); Imetal. (2002); Renzini (2003); Gao et al. (2003); 
[oiazebrook et al. (2004)), high-redshift star formation efficiency 
needs to be much higher (about 10%) in massive galaxies than in 
less massive ones. This high star formation activity at early times 
means the build-up of stellar mass does not follow the hierarchical 
build-up of total mass. Stars could form in smaller objects (not in 
the main progenitor only) with lower efficiency, provided that the 
galaxy formation process has some way to synchronise star forma- 
tion in disparate pieces. The reason being that the above mentioned 
galaxies at 2 ^ 1 — 2 have stellar populations with almost no spread 
in their stellar ages and derived star formation histories consistent 
with very short times (< 0.1 Gyr, see Macarthy et al 2004). 

(iii) There is an indication that major mergers are not the princi- 
pal drivers of star formation. 

(iv) We propose that a threshold for star formation for galaxies 
with masses ~ 10^" Mq can explain the findings above. The exis- 
tence of s uch a threshold at low red shift is well documented in the 
literature ( Martin & Kennicutt '2001'), and can be linked to feedback 
efficiency. Feedback, i.e. suppression of star formation activity, is 
expected to be very inefficient in massive galaxies, while efficient 
in less massive galaxies. 



2 BARYON CONVERSION TREE FROM SDSS 

The large size of the spectroscopic sample of the SDSS DRl pro- 
vides the means to analyse statistically properties of galaxies in- 
ferred from their spectra. Further, SDSS spectra have two important 
characteristics: a large wavelength coverage (~ 3500 A - 8500 A ) 
and a relatively high signal-to- noise (S/N 6 to 10 per resolution el- 
ement of 2A ). Previous work feauffmann et al.*2003^. ^nter et alj 
[2003; Padmanabhan et al. 2003; Kauffmann et al. 2004) accurately 
determined the total stellar mass of SDSS galaxies at the galaxy 
observed redshift. 

The observ ed spectrum of a g alaxy contains (in principle) the 
"fossil record" ^Pa nter et alj|2003l) of the galaxy's star formation 
history, that is information on the stellar mass as a function of red- 
shift; but no work has so far been done to determine the stellar mass 
formation history (the BCT). Here we attempt for the first time to 
do this. 

One important limitation is imposed by the fact that the old 
stellar populations are orders of magnitude dimmer than young 
ones, thus they have a sub-dominant effect in the galaxy spectrum. 
As a consequence, it is not possible to determine, from the observed 
spectrum, the BCT with arbitrary time resolution, especially at high 
redshift. On the other hand, we find that one can accurately deter- 
mine the BCT for time bins that are logarithmically spaced in look- 
back time. 



galaxies for this study^ . Full details of the SDSS are available at 
[http : / /www . sdss ■ orgl The spectra are top-hat smoothed to 
20A resolution for comparison with the stellar population models 
of Jimenez et al. (2003) and emission-line regions are removed ^ 
since these have a complicated dependence on the geometry of the 
ionising region and do not carry much information about the un- 
derlying stellar population. The principal strength of MOPED is 
that instead of depending on a few, possible contaminated, lines in- 
dices, it uses the whole of the rest of the spectrum in an optimally- 
weighted way, which extracts essentially all of the star formation 
history information. 

We recover the mass of stars created in 10 time bins (5Af* {ti) 
where i = 1, .., 10), which are equally spaced logarithmically in 
look-back time, separated by factors of 2.07. Table Q shows the 
centre and boundaries of the bins both in look-back time and red- 
shift. For each time bin we also recover the average metallicity of 
the stars (Z*(ti)). The final parameter is an overall dust parame- 
ter for each galaxy at the observed redshift (-Dz„); we assume an 
extinction curve as the one determined for the LMC (iGordon et alJ 
[200 J). However we are not too sensitive to the particular extinction 
curve as we have experimented with a variety of dust screens and 
found little variation in the shape of the recovered star formation. 
We use a Salpeter initial mass function, with a power-law exponent 
of X = —1.35. To model the galaxy spectra as a function of these 
param eters, we use the stellar population models'' of Jime nez et alJ 
( 2003). We should point out here that, although in principle a max- 
imum likelihood analysis could recover the star formation history 
of individual galaxies, the parameters for a single galaxy are not 
tightly constrained. In practice we statistically recover the average 
stellar assembly history from the full DRl. To do this, galaxies are 
weighted inversely with l^max, the volume in which they could be 
found and still satisfy the selection criteria for the survey. For this, 
the evolution of the stellar population and spectrum are included, 
but no size evolution is assumed. Concerns have been raised that 
the calibration of the spectra is done by using photometry with a 
larger aperture; this issue has been addressed by Glazebrook et al. 
(2003), who found that, on average, the colours from the fibres and 
from the photometry are consistent (on average), so this should not 
be a concern for the sample as a whole. To support this further, we 
show in fig.Qthe fraction of stellar mass contributed by the oldest 
three bins, for galaxies selected to be in a fixed mass range, as a 
function of the redshift of the galaxy. Encouragingly, there is no 
evidence of any significant trend towards low-redshift, where the 
aperture effects might be expected to be a problem. 

In addition to this test, we have performed fu rther checks 
on the MOPED technique. In iHeavens et alj j^OO^), we showed 
(in the supplementary information on the Nature web-site) that 
MOPED could recover the star formation correctly, given an in- 
put SFR which matched the SFR we claimed. To test this more 
thoroughly, we have generated synthetic spectra for 500 galaxies 
which have a SFR peaking at 2 = 1, corresponding to the broad 
conclusions of SFR studies pre jHeavens et alj i2004h . We also in- 
clude wavelength-dependent relative noise, characteristic of typical 
SDSS galaxies, and a systematic calibration offset at the blue end. 



2.1 Method 

The SDSS DRl main sample has apparent magnitude limits 
15 ^ ran ^ 17.77 and covers a redshift range 0.005 < 
z < 0.34, with a median redshift of 0.1. We place an addi- 
tional cut on surface brightnes s of ur < 23.0 to avoid spurious 
background contamination (see lShen et alji2003h ). leaving 96,545 



^ The calibration of the continuum blue wards of 4500 A has been im- 
proved in DR2. We will check the effect of this by analysing the DR2 sam- 
ple in a forthcoming paper. 

^ i.e., ignored in the likelihood fitting procedure. In particular the regions 
excluded are: 3700-3760 A 4840-5200 A and 6500-6800 A. 

Available at |http : / / www . physics . upenn . edu/ '^raul j| 
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Figure 1. Recovered average stellar mass fractions in the oldest three bins, 
for galaxies with masses in the range 1 — 3 X 10^^ A/q. Notice that there is 
no significant trend at low redshift, which might indicate a problem arising 
from differences between the small fibre aperture and the larger aperture 
used for the photometry. 
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Figure 3. Recovered total mass in stars as a function of redshift for the 
whole sample. 



1 .0000 r 



0.1000 ; 



0.0100 ; 



0.0010 r 



A 



0.0001 K> o. 



^ Input 

r-. Recovered 



0.01 



0.10 



1.00 



10.00 



Figure 2. Input (diamonds) and average recovered (triangles) SFR for a 
sample of 500 synthetic galaxies with wavelength-dependent relative noise 
properties typical of the SDSS. a calibration uncertainty (15% at the blue 
end, tapering to zero at 4000 Angstroms), and a variable amount of line- 
filling of the HS line (mean 0.5 of the model depth, rms 0.2). 



We also allow the HcS line to be randomly partially filled with emis- 
sion. This should be a rather thorough test that the results are not 
biased by noise, calibration or line-filling. The input and average 
recovered star formation rates are shown in Fig. |2| which illustrates 
that indeed MOPED does recover the input star formation history 
without any biases. 

Although we are mostly interested in the 10 5M^.(ti), in total 
we have 21 parameters which we want to constrain from 96,545 
SDSS spectra (each of which has about 2000 elements) using a 
maximum likelihood approach. This is an extremely computation- 
ally expensive task, which we can handle by resorting to 1) a data 
compression algorithm and 2) the now widely used Markov Chain 
Monte Carlo (MCMC) technique. 

The data compressio n algorithm MOPED 

iHeavens. Jimenez & Lahavl l2000l) enables us to explore effi- 
ciently the parameter space and place error bars on the recovered 
parameters. The algorithm linearly combines the 2000 flux 
elements in each spectrum in 21 MOPED coefficients, one for 



each parameter, that contain all the relevant information (i.e. 
the method is lossless). F or further details about MOPED see 
iHeavens. Jimenez & Laha\it2000) . We have already demonstrated 
the usefulness of MOPED and the MCMC implementation 
to recover galaxy physical properties from different galaxy 
spectroscopic samples, includ ing the SDSS EDR and DRl 
i Reichardt, Jimenez & Heavens 2001: Panter. H eavens & Jimenez! 
t2003.,: Heavens. Panter. Jimenez & Dunlon ,2004h . Essentially, it 
works by weighting the flux data in a way which preserves the 
information inherent in the complete spectrum, so basically no 
information is lost and the procedure is almost equivalent to using 
all of the flux values in the spectrum. For an individual galaxy, 
there may still be residual degeneracies in the solution, but we 
have demonstrated with tests that for a very large sample, such as 
in SDSS, the average solution recovers the input extremely well. 

The SDSS survey is a magnitude-limited survey and here we 
aim at deriving average properties, i.e. the stellar assembly his- 
tory, of galaxies as a function of their present-day stellar mass. In 
principle, this could introduce a bias since for a given mass, galax- 
ies with young stellar population, tend to be brighter than galaxies 
with an old stellar population; thus blue galaxies which form stars 
recently, would be preferred over passively evolving ones. However 
this does not matter for our purposes for the following reasons: a) 
galaxies have been effectively selected by mass from a complete 
redshift and magnitude-limited sample. This can be seen in Fig.|3| 
where the recovered total mass is plotted as a function of the ob- 
served redshift. Note that as expected, for a mass-selected sample, 
less massive galaxies are at low redshifts while most massive galax- 
ies occupy the whole redshift range; b) we will determine galaxy 
properties averaged in mass-bins; cj the recovered mass fractions 
for each time bin for each galaxy are properly shifted according to 
the observed redshift of the galaxy and assuming that the size of 
galaxies does not evolve over the period of time covered by the ob- 
served redshift range. Note that the results on stellar mass fractions 
are insensitive to the redshift of the galaxies used - see Fig.|3l To 
determine average properties as a function of mass we average our 
results in mass with mass bins equally-spaced logarithmically by 
factors of ten. 
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Table 1. Centre and boundaries of the bins in redshift (z) and look-back time tiookback in Gyr 



bin# lower 2 centi'e 2; upper 2 lower tiookback(Gyr) centre tiookback (Gyr) upper tiookback(Gyr) 



1 


0.0007 


0.001 


0.00145 


0.00966 


0.014 


0.0200 


2 


0.00145 


0.0021 


0.003 


0.0200 


0.029 


0.0414 


3 


0.003 


0.006 


0.0063 


0.0414 


0.06 


0.0857 


4 


0.0063 


0.012 


0.013 


0.0857 


0.12 


0.1776 


5 


0.013 


0.0179 


0.027 


0.1776 


0.26 


0.3677 


6 


0.027 


0.0419 


0.057 


0.3677 


0.53 


0.7614 


7 


0.057 


0.0839 


0.125 


0.7614 


1.10 


1.5767 


g 


0.125 


0.186 


0.287 


1.5767 


2.27 


3.2650 


9 


0.287 


0.456 


0.786 


3.2650 


4.70 


6.7609 


10 


0.786 


1.755 


5.000 


6.7609 


9.7 


12.000 




Figure 4. Top panel: fraction of stellar mass of the galaxy formed in dif- 
ferent time bins as a function of total stellar mass. Dashed line corresponds 
to the oldest bin (9.7 Gyr), dotted line to the second oldest (4.7 Gyr) and 
solid line to the third oldest (2.3 Gyr). The other solid lines correspond to 
younger time bins (1.1, 0.53, 0.26). Bottom panel: fraction of the observed 
total stellar mass of the galaxy created as a function of time for different 
galaxy stellar masses (determined at the observed redshift). 



Figure 5. Top panel: fraction of the final dark matter mass that has been 
virialised in the most massive progenitor, as a function of the total stellar 
mass of the galaxy (at the observed redshift). The different lines con'espond 
to the amount of dark matter that is virialised in the given time bin. Bot- 
tom panel: fraction of the final dark matter mass that has been virialised in 
the most massive progenitor as a function of tiookback ■ The different lines 
con'espond to different total stellar masses of the galaxy at the observed 
redshift. 



2.2 Results 

We find that more massive galaxies transform their gas into stars 
earlier than less massive ones. The top panel of Fig. |4] shows for 
each time bin how many baryons (as a fraction of the observed total 
stellar mass) were converted into stars, fs , as a function of observed 
stellar mass — M,{t = 0). The different lines correspond to 
different time bins: dash oldest, dots second-oldest and continuous 
lines denote younger bins. The bottom panel shows the stellar mass 
assembly history 5Mt{ti) for different observed stellar masses. 



There is a clear correlation between baryon conversion effi- 
ciency and present-day stellar mass. This can be see from Fig.|4] 
In galaxies with M* larger than 3 x lO^'^ M©, more than 60% of 
their present stellar mass was already in place at redshift 1.7. In 
terms of look-back time it means that the stars of massive galaxies 
(i\/* ~ 10^^ Mq) were essentially formed (if not in place) more 
than 9 Gyr ago, or just ~ 4 Gyr after the big bang. Conversely, 
for _A/« < 10^ M0, more than half of their stellar mass remains 
unformed at 2; = 0.2 or ~ 3 Gyr ago. 
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3 DARK MATTER AND BARYON ASSEMBLY HISTORY 



3.1 The number of progenitors of galaxies 



We now compare the stellar formation history we just recovered, 
with the dark matter assembly history. To compute the mass of 
dark matter halos as a function of time we use two approaches: 
first, we generate multiple realizations (10^) of the merger history 
of dark matter halos of several masses in the ra nge 10^ to 10^'^ M^ti. 



This i s done using the algorithm described in lSomerville & KolattI 
j 19991) for t he standard LCDM co smology (fio = 0.27, Ao = 0.73, 
h = 0.71, ISoergel et alj i2003h ). This algorithm reproduces the 
merger histories of halos seen in N-body simulations of structure 
formation I Somerville et al. 2000; Wechsler et al. 2002), especially 
at low redshift, which is the range of interest. One free parameter 
in the algorithm is the value of the mass that is considered accreted 
instead of merged. Agreement with CDM simulations is achieved 
when everything below 1% of the final halo mass is considered 
accreted and this is the value we use. Secon d, we use the fitting for - 
mula for the mass accretion history from IWechsler et al] i2002h . 
This is obtained fro m numerical N-body simulations performed 
with the ART code iKravtsov et alJll997h . We find that both ap- 
proaches show the same qualitative behaviour illustrated in Fig.|5] 
howe ver since the Extended Press -Schechter approach, at the base 
of the lSomerville & KolattI i 19991) algorithm, is not a perfect fit to 
N-body simulations, especially at high-redshift and for very mas- 
sive halos, we will present here only results obtained using the sec- 
ond approach. 

We compute the dark matter mass of the most massive progen- 
itor that is virialised in each redshift bin as a function of the total 
mass of the dark halo. The top panel of Fig. |5| shows the fraction 
of the final dark matter mass that has been virialised in the most 
massive progenitor, /d , in a given redshift bin (line styles same as 
in Fig. |4j as a function of the final stellar mass in the dark halo. 
The stellar mass of the dark halo is obtained from the dark one as- 
suming the universal baryon fra ction as de t ermin ed by WMAP 
(ncDMM = 4.8 = /DM/b;|S2i^eLi] 1200?)), and that at 
2^0 only about 6% of the baryons are in stars iFukugita 2003). 

The bottom panel of Fig. |5| shows the mass assembly history 
of the dark halo. 

A comparison of Fig. |4|and Fig. |5] indicates that dark mat- 
ter assembly and formation of stars do not follow each other. For 
example, for A'/« > 10^^ Mq less than 50% of the dark matter 
is assembled in the main progenitor at 2 > 1.7 (iiookback ~ 9.7 
Gyr), while more than 75% of the stellar mass is already formed. 
On the other hand, for stellar masses smaller than lO" Mq, 20% of 
the stellar mass is formed in the same time bin while already 60% 
of the dark matter is in place. This hints at a role of early stellar 
feedback in these halos as we discuss later in §4. 

While in the hierarchical LCDM model for structure forma- 
tion the more massive CDM structures form late, we find obser- 
vational evidence for early star formation of giant galaxies. This 
can happen because of two reasons: a) massive galaxies are formed 
by mergers with smaller ones, each carrying an evolved stellar 
population b) these massive galaxies are already in place at high- 
redshift and the dark matter halo has been assembled at the same 
time as the stellar population. If the dark matter halo collapse trig- 
gers star formation then one would expect a) to be the case; how- 
ever observations of e.g ., old elliptical galaxies at high redshif t 
(e.g.lDunlop et alJJl996l) ; lNolan et alj ^^-^^^^ i2004) ; 
ISaracroe^ni2003l) ) seem to support the second scenario, at least 
in some cases. 



It is clear that since massive dark halos are predicted, on average, 
to assemble later than the stellar population they contain, the stars 
may naturally have formed previously in smaller dark halos that 
subsequently merged. By considering the (dark) mass accretion his- 
tory of the most massive progenitor of a halo, and the star-formation 
rate, we can constrain the minimum number of progenitors form- 
ing a halo, and the star formation efficiency as outlined below. In 
other words, if the most massive progenitor at a given redshift z 
carries enough baryons to form all the stars that should be in place 
by z then only one progenitor (the most massive one) is needed. If 
the minimum number of progenitors forming a halo is larger than 
1, massive galaxies must have form by mergers with smaller ones, 
each carrying an evolved stellar population. However, if the min- 
imum number of progenitors is 1, then the LCDM paradigm can 
accommodate old elliptical galaxies at high redshift. 

For example, consider a galaxy with stellar mass Af* — 
M0 (see Fig.|4]and Fig.|5}. In the oldest time bin 76% of the stellar 
mass is already in place (Af, (tio) = 0.76A/,), yet the virialised 
fraction of the dark matter halo is less than 50%. If we assume case 
a) (that is star formation happens uniformly in all halos and sub- 
halos that ultimately end up forming the final galaxy, and that 100% 
of the dark matter is virialised at all times) then the fraction of the 
total baryonic mass converted into stars was 4.6%. On the other 
hand, observations of old elliptical galaxies at high redshift can be 
explained within the LCDM paradigm if we assume that the main 
virialised progenitor harboured all the stars observed, in this case 
the fraction of mass converted to stars in the dark matter halo must 
have been ~ 10%. This requires only a modest enhancement in 
star formation efficiency at high redshift, but close to the maximum 
efficiency observed in giant molecula r clouds today, which is < 
10% (e.g. JPadoan & Nordlunjiiool ). 

More specifically, we can assume that the stellar mass 
Af, (tio), was in progenitors (the most massive progenitor and pos- 
sibly other sub halos), whose cumulative dark matter must have 
been at least A/,(ti()) x foM/b/f, where / ^ 1 parameterises 
the star formation efficiency, and the fraction of baryons that gets 
turned into stars, and depends on the mass of the sub-halos. We 
approximate /(Afdark) as /s(A^dark//DM/b x 0.1) (as 0.1 is the 
minimum efficiency in the oldest time bin). 

The minimal number of progenitors can thus be obtained by 
minimisation as a function of the sub halo dark mass Msh'- 



N{Msii) = 1+max 



.fDM/bM.(tlo)//(AfsH) - Afdark(tio) 

Msh 



We obtain a minimum number of 1 and / ~ 10%, in agreement 
with the more heuristic argument above. If the main progenitor did 
not have enough baryons to accommodate A/»(tio), then the mini- 
mum number of progenitors would be greater than one. If the mini- 
mum number of progenitors is one and the star formation efficiency 
is constrained to be reasonable, then all the old stellar population 
could have been formed in the main progenitor. 

Thus this suggests that in the LCDM paradigm, the massive 
old galaxies that we see today could have been made of merg- 
ers of few progenitors, each carrying an old stellar population, 
in agreement with other indications that elliptical galaxies are al- 
ready formed_at_z_ > 1 (e.g. Bower et al. (1992); Peacock et alj 



199S) 
2002) 



Lilly et alJ <1998l); iBrinchmann & Ellisi (.2000.); .Im et all 
Renzinil |2003l) ; iGao et alj llOO^V . I Saracco et al 1 <2003h ; 



iGlazebrook et a l] j2004 VA small number of mergers can also nat- 
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Figure 6. Ratio of mass transformed into stars in a galaxy to the baryons 
that are accreted onto the main progenitor, which typically contains > half 
the mass. Crosses correspond to a stellar mass of 10^^, asterisks to 10^^, 
diamonds to 10^", triangles to 10® and squares to 10* Mq. A value for 
a!Af*(t)/iiJ\/baryons(t) Smaller than 1 means that not all the baryons (the 
nucleosynthesis value as recently constrained by WMAP) have been trans- 
formed into stars, while a value larger than 1 indicates that gas previously 
available in the galaxy has been turned into stars: either the recent accre- 
tion has triggered star formation in the main progenitor, or it is going on 
elsewhere. Note that for stellar masses below 10^" Mq, the star formation 
efficiency peaks at iiookback ^ 

1-2 Gyr 



urall y explain the tightn ess of the observed colour-magnitude rela- 
tion jBoweretalll992h . 



4 TIME EVOLUTION OF STAR FORMATION 
EFFICIENCY 

For each of the time bins we compute the ratio of the newly-formed 
stellar mass to the baryonic mass added to the main progenitor, 
assuming the nucleosynthesis baryon fraction. Fig. |6| shows the 
above ratio as a function of look-back time (tiookback) for differ- 
ent masses: crosses correspond to a stellar mass of 10^^ Mq, as- 
terisks to 10" M©, diamonds to lO"' M©, triangles to lO'' M© 
and squares to 10* M©. This is a measurement of how much gas is 
transformed into stars as a function of the newly-added baryons. A 
value of 1 indicates that the mass of baryons accreted to the main 
progenitor matches the mass converted into stars. A value higher 
than 1 shows that the accretion or merger was accompanied by a 
greater mass of triggered star formation somewhere in the galaxy. 
Fig. |6| clearly shows that for stellar masses above 10^" M© this 
ratio is never greater than 1. Clearly we are comparing mass ac- 
creted on to the main progenitor with stars created in any of the 
progenitors, and this should be borne in mind in interpreting Fig.|5| 
However, since the main progenitor contains ~ 50% of the final 
mass even at a lookback time of 10 Gyr, and this fraction is weakly 
dependent on mass, the efficiency of conversion of baryons to stars 
in the galaxy as a whole is unlikely to alter the main conclusions of 
Fig.|5| where the differences between objects of different present- 
day masses typically far exceed a factor of 2. 

For the most massive galaxies at early times this measure of 
star formation efficiency is close to 40%. For stellar masses below 



10^"M©, the efficiency is only of about 6-8% at tiookback = 10 
Gyr, but grows to 100% at tiookback = 2 Gyr, and then decreases. 

For galaxies with stellar mass smaller than 10^ M©, this in- 
crease in star formation efficiency rises until fiookback ~ 0.5 Gyr, 
at which point it reaches 300% efficiency, which means that more 
gas is transformed into stars than the baryons brought into the par- 
ent dark halo by accretion. This points toward a picture where these 
low-mass gas-rich galaxies see a lot of their gas reservoir trans- 
formed into stars due, for example, to a merger or accretion event. 
Another possibility is that star formation is proceeding rapidly in 
other sub-halos, which subsequently merge with the parent. This 
scenario is not supported by the merger histories of low mass halos 
(e.g. Sheth &l Tormen 2004; Wechsler et al. 2002) since small halos 
at present time have almost no merging from smaller subhalos. 

Thus massive galaxies have a high star-formation efficiency 
at early times and then evolve "passively", with fresh infall of gas 
being suppressed or turned into stars with low efficiency, possibly 
because it is likely to be too hot. Small galaxies seem to accrete 
mass passively at early times and form stars very efficiently later. 

Conversely, the probability distribution for dark halo merger 
events peaks at higher redshift for small halos and at lower redshift 
for large halos. In a A-dominated Universe, merger probability is 
suppressed at 2: < 1 especially in low density regions, where small 
galaxies are most likely to be. Thus there seems to be no correlation 
between halo virialization or dark matter merger events and star 
formation efficiency. 

However, one could imagine explaining this trend, for exam- 
ple, by postulating the existence of a "threshold" for star forma- 
tion: once this threshold is crossed, all available baryons are turned 
into stars (as in an "infall model"), then afterward galaxies evolve 
passively (as in a "closed box model"). In this toy model, if this 
threshold is crossed at very early times in the progenitors of mas- 
sive galaxies, one would expect these galaxies to form stars very 
efficiently early on then evolve passively. For progressively smaller 
galaxies this threshold is met at increasingly later times. 

An alternative explanation is that stellar feedback is the re- 
sponsible for the lack of star formation in small galaxies at early 
times. Since the escape velocity in these systems is smaller than 
in more massive galaxies, gas can leave the dark matter halo more 
easily. Only the very massive systems are able to retain their gas 
and convert a majority of their gas into stars. If so we find that 
M* ~ lO^'^ M© is the characteristic mass that defines the border 
between efficient and inefficient feedback. 



5 CONCLUSIONS 

We have determined for the first time the baryonic conversion tree 
for galaxies. We were able to do this with observations at 2: < 0.3, 
using 96545 galaxies from the SDSS DRl spectroscopic sample. 

In the hierarchical structure formation model, massive dark 
halos form (i.e. virialize) later than smaller halos, from mergers of 
smaller units (e.g. Lacev & Cole ( 1993, 1994); Lin et al. (2003)). 
This model has been thoroughly tested against numerical cold dark 
matter simulations. Naively one could expect that the dark matter 
halo collapse should trigger baryonic gas transformation into stars; 
in addition subsequent dark matter mergers should produce star for- 
mation episodes. 

Instead, for the the stellar assembly history we find that: the 
more massive galaxies have old stellar population and massive, old 
elliptical galaxies are already in place at 2: ~ 1. This has been 
known for a long time and sometimes it is referred to as "down- 
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sizing , JCowie etalJll996l) . We find that massive galaxies have 
transformed more gas into s tars at higher redshif t (in agreement 
with high z observations e.g.. lKodama et all i2004) l and then star 
formation was suppressed, while less massive galaxies transform 
more gas into stars at low redshift. 

So one should not be surprised to see abundant red ob- 
jects at high redshift in the LCDM paradigm, these objects can 
form in virialised h a los if star formation efficiency is high. In- 
deed, I Jimenez et alJ i 19991) have shown that single-halo hydro- 
dynamical models would require an increased star formation ef- 
ficiency for more massive galaxies, higher than the few per- 
cent found today in giant molecular clouds, in agreement with 
the value determined from the "fossil record" in the present 
work. Our findings, based only on observations at z < 0.35 
(the "fossil record"), are in agreement with a suite of inde- 
pendent, z > 1, observations: observations of old elliptical 
galaxies at high r edshift (DunloD et al. 1996; SDinrad et al. 1 997t 
iNolan et aljboOSh . ind ications that elliptical galaxies are alread y 
formed at z > 1 



(e.g. 'Bower et alj il993): iPeacock et alJ jl998[) : 

i Lillv et alJ il998t): Is rinchma nn & Ellis! jaOOOl): llm et alj 120021) : 
kenzini' 120031): loao et all l2003h : iGlazebrook et alj J2004) ). It 
also nicely explain s the tightness of the colour-magnitude relation 
jBower et aljl992h . 

On the other hand, we find that small galaxies seem to accrete 
mass passively at early times and see a lot of their gas reservoir 
transformed into stars at late times. Since the probability distribu- 
tion for dark halo mergers peaks at low redshift for massive ha- 
los and high redshift for small halos, we conclude that dark mat- 
ter mergers and star formation are not correlated. We speculate that 
one possibility to explain the apparent and illusory anti-hierarchical 
nature of the stellar assembly history is the existence of a thresh- 
old for star formation: once the threshold is crossed all available 
baryons are turned into stars ("infall model") and afterward galax- 
ies approximatively evolve passively. The threshold is met at very 
early time fo r massive galaxies and a la ter time for less massive 
ones (see e.g. lHeavens & Jimen ez 1 1999)). 

A star formation thre s hold h as been observed in disk galax- 
ies bv lMartin & Kennicutd i200ll) and there has been some recent 
additional evidence from the formation of dust lanes in disk galax- 
ies iDalcanton et al. 2004) that this threshold may t ake place at 
Ve ^ 100 km s~^, in agreement w ith the findings of IVerde et alJ 
j2002h and lKannappan et alj liooj) who also found a transition at 
about 100 km s~^ for star formation efficiency. This Vc value cor- 
responds to the characteristic mass found here (Af. ~ lO" Mq), 
that defines the border between efficient and inefficient star for- 
mation. This characteristic mass has been related to feedback effi- 
ciency threshold (e.g., Dekel & Silk 1986, Dekel & Woo 2003 and 
references therein). 

As we do not yet have a fundamental theory for galaxy 
formation and given the complexity involved in studying the 
process with hydrodynamic-N body simulations, we hope that this 
new determination of the baryonic conversion history will be a 
useful observable to gauge galaxy formation models against. 
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